[hep-ph/0 104246 

CERN/2001-114 
Edinburgh 2001/02 
RM3-TH 01/5 



Extrapolating Structure Functions to Very Small x 

Guido Altarelli/ Richard D. Ball^' ^ and Stefano Forte^-g 

^Theory Division, CERN 
CH-1211 Geneva 23, Switzerland 

^Department of Physics and Astronomy, University of Edinburgh, 
Mayfield Road, Edinburgh EH9 3JZ, Scotland 

^INFN, Sezione di Roma Tre 
Via della Vasca Navale 84, 1-00146 Rome, Italy 



Abstract 

We review small x contributions to perturbative evolution equations for parton 
distributions, and their resummation. We emphasize in particular the resummation 
technique recently developed in order to deal with the apparent instability of naive 
small X evolution kernels and understand the empirical sucess of fixed-order per- 
turbation theory. We give predictions for the gluon distribution and the structure 
functions F2{x,Q'^) and Fi{x,Q'^) in an extended kinematic region, such as would 
be relevant for THERA or LEP+LHC ep colhders. 
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Measurements of the inclusive structure functions F2{x,Q^) and Fl{x,Q^) at HERA 
have shown that the scahng violations of structure functions are in extremely good agree- 
ment with the perturbative next-to-leading order (NLO) QCD prediction, down to the 
smallest values of x, and for all ^ 1 GcV^ [1]. This agreement is surprising in that 
it is known that perturbative corrections beyond NLO in are enhanced by powers of 
^ = ln(l/a;), and thus one would expect higher order corrections to be sizable whenever 
c^siQ'^)^ ~ 1; i-e. in most of the HERA kinematic region. Whereas techniques for the 
inclusion of small x contributions to leading twist evolution equations have been known 
for some time [2,3], only recently did a consistent picture of the general structure of these 
contributions and their resummation emerge. Indeed, considerable theoretical progress 
has been spurred by the determination [4] of next-to-leading corrections to the BFKL ker- 
nel, which allows the computation of the next-to-leading log(l/a;) (NLLx) contributions 
to anomalous dimensions to all orders in ag. Specifically, it is now understood that the 
inclusion of NLLx contributions leads to instability [5] of perturbative evolution, unless it 
is suitably combined with a resummation of the collinear singularities [6-8] which are re- 
summed order by order in the standard QCD evolution equations. Furthermore, the NLLx 
perturbative corrections give rise to increasingly large contributions to high orders of per- 
turbation theory [9,10] that make a nonsense of the perturbative expansion and call for 
an all-order resummation of the small-x behaviour of the anomalous dimensions [11,12]. 

Practical methods to deal with these issues have been developed recently [7, 13], and 
lead to a resummation prescription which is amenable to numerical treatment and direct 
comparison with the data. It then appears that the observed smallness of perturbative 
higher order corrections at small x can be accommodated within the current knowledge of 
the general structure of anomalous dimensions, but it poses very stringent constraints on 
the form of the unknown higher order terms. Furthermore, even when these constraints 
are respected, so that, as required by the data, deviations of the behaviour of the observ- 
able structure functions from the fixed next-to-leading order prediction are very small, 
still non-negligible modifications of the fitted parton distributions at small x are found. 
This, because of ambiguities in the resummation procedure, entails larger uncertainties 
on parton distributions at small x. Likewise, these corrections have a sizable impact on 
the extraction of from small x data, both on the central value and the estimates of 
overall theoretical uncertainties [13,14]. 

In the wider kinematic region available at THERA the small differences between re- 
summed and fixed-order predictions could be put to more stringent tests. This would 
allow one to pin down more precisely the ambiguities in the resummation procedure, 
thereby reducing the uncertainty on parton distributions at small x and on precision de- 
terminations of as at small x. Also, the possibility of reaching smaller values of x for given 

would allow a test of resummed perturbation theory in a region where the relevant re- 
summation parameter is large, and also to see whether the perturbative description of 
scaling violations remains satisfactory or starts to break down, as is often suggested [15]. 

Here we briefly review our current understanding of resummed perturbation theory at 
small X. We then give predictions for the gluon distribution and the structure functions 
F2 and Fl in two different resummation scenarios, and compare these to fixed next- 
to-leading order results in the kinematic range which is relevant for THERA. This is 
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essentially the same kinematic region accessible at a hypothetical lepton-hadron collider 
obtained combining LEP with the LHC, so our predictions would also be relevant at such 
a machine. 



1 Duality of small x evolution 

The basic result which allows the determination of contributions to anomalous dimensions 
which are logarithmically enhanced in x to all orders in the coupling, and thus their 
inclusion in evolution equations, is the duality of perturbative evolution [7,16]: because 
leading-twist evolution of structure functions takes place both in x and Q^, it admits 
a dual description in terms of equations for evolution in t = ln(Q^/yU^) or evolution in 
^ = ln(l/a:;). This property is easy to prove [16] when the coupling is fixed, and can be 
shown to remain valid when the coupling runs by explicit order-by-order perturbative 
computation [12]. 

Let us first consider for simplicity the case (relevant in the very small x limit) of a 
single parton distribution G{^, t), identified with the dominant eigenvector of perturbative 
evolution. The pair of dual evolution equations are then 

jG{u) = p(e,«s)®G(e,t), (1) 

= ir(t,a.)®G(e,t), (2) 

The convolutions on the right-hand sides of the dual evolution equations are with 
respect to ^ in the first equation as) is the usual splitting function) and with respect 
to t in the second equation; as = as{t) and is unaffected by convolutions. Duality means 
that the solutions to these equations coincide up to higher twist corrections provided the 
respective boundary conditions and kernels are suitably matched. 

The detailed form of the matching of boundary conditions is irrelevant for our purposes, 
but it is important to notice that the matching is such that the boundary condition to (|I|) 
depends only on ^ (and not on t) and the boundary condition to (^ depends only on t 
(and not on ^) as required by factorization. The matching of the kernels is given by the 
duality equation 

x(7(iV,«s),«s) = iV, (3) 

or equivalently its inverse 

7(x(M,«,),a,) = M. (4) 

Here 7 is the usual anomalous dimension, related to the splitting function by Mellin 
transformation with respect to ^: 

POO 

7(iV,«,)= / d^e-^^ P{^,as). (5) 
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The relation between x{M, as) and K{t, as) is somewhat more comphcated because, upon 
Melhn transformation with respect to t the running couphng as{t) on the right-hand side 
of Eq. (|l]) becomes a differential operator. The relation between the evolution kernel 
K{t, as) and the dual kernel xi^y cts) can nevertheless be determined order by order in 
perturbation theory [12]: defining 

K{t,as) = asKoit) + alK,{t) + . . . , (6) 
x{M,as) = asXo{M) + alxi{M) + . . . , (7) 

we get 

oo 

Mt 



Xo(M) = / dte-^^' Ko{t); 

J — oo 

" ' i-co 4,r2 y:'„'(M) ^ ' 

where Pq = — |n/ is the first coefficient of the QCD P function. 

It follows from the form of the duality equation @ that knowledge of the leading 
(next-to-leading, . . . ) term in the expansion of x in powers of as at fixed M determines 
the leading (next-to-leading, ... ) term in the expansion of 7 in powers of at fixed 
as/N: i.e. defining further 

7(iV,a,) = 7,(^)+a,7,,(f) + ... , (9) 



then 

xo(7s(^)) = -, %sm = - ^:)'')z:: , ■■■ ■ (lo) 



Likewise, knowledge of the leading, next-to-leading, . . . terms in the expansion of 7 in 
powers of as at fixed determines the leading, next-to-leading, . . . terms in the expansion 
of X in powers of as at fixed ag/M: writing 

7(iV,«,) = «,7o(iV)+«hi(^) + ---, (11) 
xiM,as) = X.(t) + «.X..(t) + --- , (12) 



then 

7o(x.(t)) = -, xss(t) = - :;r;f;; , ••• • (is) 



M ^ 7i(x.(!f)) 



It should be understood that the running coupling corrections Eq. (H) are always included 
in the definition of x in the above equations. 

Because the ^ evolution equation is essentially the same as the BFKL equation (up to 
factorization scheme and scale choices, which become relevant beyond leading order [4,17]) 
the duality relation can be viewed as a consistency condition between this equation and 
the standard renormalization group equation for moments of structure functions in the 
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Figure 1: Plots of different approximations to x- fhe BFKL leading and next-to- leading 
order functions (0), a^xo and a^xo+^sXi (dashed); the LO and NLO dual a^Xs and asXs+ 
a^gXss ( |12|) of the one and two loop anomalous dimensions (solid), and the double-leading 
functions at LO and NLO defined in Eq. (|T5|) (dotdashed). All curves are computed with 
a, = 0.2. 



region of their common validity (i.e. large and small x). Hence, knowledge of the 
BFKL kernel K{t, as) (§) can be translated into information of x (0)? which in turn can 
be used to gain information on the logarithmically enhanced contributions 7s, 7ss, . . . (^) 
to the anomalous dimension 7(05, A^), and conversely. In fact, the leading-order equation 
in ( p^ has been known for a long time [18]; the new insight here is that this is just a 
consequence of a more general duality. 



2 The Double— Leading expansion 

Only the first two orders in the expansion of x at fixed M and 7 at fixed N are currently 
known. While the perturbative expansion of 7 is well-behaved, in the sense that agji is 
a small correction to 70 for reasonable values of the coupling constant, the perturbative 
expansion of x is very poorly behaved, in that the NLO correction xi completely changes 
the qualitative shape of the kernel. In particular (see Fig. |l]), in the physical region 
< M < 1 the LO kernel has simple poles with positive residue at M = and M = 1 
and a minimum in between. The NLO correction xi instead has higher order poles with 
negative coefficient, and, for any realistic value of as (essentially, for all as ^ 0.03) the full 
NLO function has just a maximum (for smaller as it has a minimum and two maxima) [9]. 
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It is easy to show that the solution to the evolution equation determined by a kernel with 
this shape displays unphysical oscillatory behaviour in the limit as a; — >^ 0, and thus, in 
particular, leads to negative cross-sections [5]. 

Because the Mellin transform (g) of = ln*^(QV/i2) is {k - l)!/M^-\ the presence 
of 1/M poles in the kernel x is related to collinear singularities: indeed, according to 
Eq. (^) the coefficients of these singularities are determined by knowledge of the anomalous 
dimensions 70, 71, . . . in the usual renormalization group equations, which resum collinear 
singularities. It is easy to understand [7] why these singularities lead to a series of poles 
in M = with alternating signs. Indeed, recall that momentum conservation implies 
that the largest eigenvalue of the anomalous dimension matrix vanishes at iV = 1, i.e. 
7(1, as) = 0, which by duality (|^) implies x(0,as) = 1- It follows that if, in the vicinity 
of M = 0, Xs behaves as 

y ^ = ^-^ + ^+ ■ (u) 

the series of poles in Xs, Xss, ■ ■ ■ actually sums up to the regular behaviour x(0, a^) = 1. 

The poles in x as M ^ are summed to all orders into Xs, Xss, ■ ■ ■ , and thus the 
undesirable behaviour of the expansion of x can be removed by defining order by order an 
improved expansion. Namely, we define a double leading expansion where to each order 
in as both the terms present in the expansion in powers of as at fixed M and at fixed 
as/M are included: 

X(M, as) = Kxo(M) + Xs (t) - ^] 

+as Kxi(M) + Xss (t) - (ife + fe) - /o] + • ■ • • (15) 

In this expansion, in the vicinity of M = the singularities of xo, Xi, ■ ■ ■ are resummed 
into Xs, Xss, while the subtraction terms avoid double-counting of these contributions. 
Note (see Fig. |I]) that at larger values of M the shape of Xo; Xi, • • • is reproduced, but in 
most of the M range the kernel Eq. (|T5|) coincides with the (dual of) the standard anoma- 
lous dimensions 70 and 71, consistent with the empirical smallness of small-x correction to 
perturbative evolution. This also has the significant implication that the double-leading 
expansion of x is as stable as the usual expansion of 7 at fixed N. 

It is easy to show that the corresponding double leading expansion of 7, 

7(iV,a,) = [a,7o(iV)+7s(t)-^] 

+as [«s7i(iV) + 7ss (t) - «s + t) - eo] + • • • , (16) 

is consistent with duality, in that x (0) and 7 (|16D are dual to each other order by order 
in the double-leading expansion, up to higher order corrections. Hence, for practical 
applications we may directly use the double-leading anomalous dimension (|l^) in the 
usual evolution equation (|l|). This will ensure that collinear singularities are resummed 
according to the renormalization group in the usual way, while leading logs of 1/a; are 
consistently included up to next-to-leading order. 

For actual phenomenology, the full set of anomalous dimensions and coefficient func- 
tions are needed. It is easy to see that the double-leading expansion is consistent with 
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diagonalization of the anomalous dimension matrix, in the sense that one may equiv- 
alently, up to subleading corrections, construct a two by two matrix of double-leading 
anomalous dimensions and diagonalize it, or else construct directly a double-leading ex- 
pansion of eigenvalues and projectors. Because one of the two eigenvectors of 7 is free 
of small-x singularities, so its double-leading expansion coincides with the standard ex- 
pansion at fixed N, the latter procedure is in practice simpler. Hence, the double leading 
expansion can be fully defined in terms of the expansion of the large anomalous dimension 
eigenvalue, and of the quark-sector matrix elements which determine the projectors on 
the eigenvectors. Likewise, one can construct double-leading coefficient functions, and 
prove that the expansion transforms consistently upon changes of factorization scheme. 
Detailed proofs and results needed for a practical implementation are given in ref. [13]. 



3 Resummation 

Even though the difference between the double-leading expansions of x dUD and 7 ([l6|) 
is subleading, it can in practice be large when M ~ 0.25. Indeed, recall that duality (^) 
implies x = It is clear from Fig. |l| that in the region M ~ 0.25 the difference between 
the leading order and next-to-leading order double-leading curves is small for any fixed 
value of M, but it is quite large for a fixed value of x = N, because the curves are almost 
parallel to the M-axis: the LO BFKL curve has a minimum at M = 1/2. Since 7 is a 
function of N, in this region the perturbative solution (|10|) of the duality relation is 
not good and the expansion of 7 is not well behaved. 

A possible way out is to determine the double-leading 7 (0) from the double-leading 
X ( |T5| ) by solving the duality relation @ exactly (rather than perturbatively). This 
can be done for instance by numerical methods, or equivalently by differentiating with 
respect to t the solution of the evolution equation (||) determined using the double- 
leading X kernel (ITH), as in ref. [8]. However, this approach, besides being cumbersome 



to implement in standard evolution codes, has the shortcoming that it hides a genuine 
perturbative ambiguity. Indeed, in this way the perturbative expansion of 7 is in practice 
stabilized by assuming that in the region M ^ 1/2 the (large) subleading corrections to 7 
will be such as to reproduce the shape of x, as computed to some fixed perturbative order, 
or possibly further improved according to a model of its behaviour at large M ~ 1 [8]. 

We instead prefer to use only the available perturbative information on 7, without 
making model-dependent assumptions. It can be shown [12] that the poor perturbative 
behavior of the expansion of 7 at fixed ag/N manifests itself in a rise of the associate 
splitting functions: Pgs/Ps ~ c^sC,, Psss/Ps ~ a^^^ and so on. This rise can be removed 

5 — i-OO ^ — >oo 

by simply subtracting at each order a suitable constant Cj from Xi (computable order 
by order in perturbation theory as a function of Xi and their derivatives at M = 1/2), 
and then determining jss... from the subtracted Xi- Thus, the expansion of 7 (|^ can be 
stabilized by just reorganizing the perturbative expansion of x'- 

x{M,as) = a,Xo{M) + alxi{M) + . . . (17) 
= asMM) + alxiiM) + . . . , (18) 
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where 



a,Xo{M,a,)=a,Xo{M)+AX, Xi{M) = Xi{M) - a, (19) 
for z = 1, 2, . . . , and thus 

oo 

AA = J]<+V (20) 

n=l 

If X has a minimum, then its value at the minimum coincides [7] with the value of xo at 
its minimum M = 1/2, namely 

A^Xo(|)=Xo(i) + AA. (21) 

Since the value of x at its minimum determines the asymptotic behaviour of the structure 
function as x ^ 0, this implies that in order to remove the perturbative instability it is 
necessary and sufficient to resum the asymptotic small x behaviour into the leading order 
kernel xo- The perturbative instability signals the fact that the all-order asymptotic 
behaviour must be known to all orders. 

Of course, we are free to use any particular truncation of A A (^Oj): for instance, we 
could simply take x to coincide with its NLO form in the double-leading expansion. 
Eq. ( p!8D then provides us with a stable perturbative expansion of 7, which at NLO is 
very close to the exact dual of x, the large subleading corrections having been resummed 
in a minimal way. In this way Eq. ([T8|) gives us a simple prescription which completely 
stabilizes the double-leading expansion of 7 whenever the double-leading expansion of x 
is also stable. Hence, any specific resummation of x (such as that constructed in ref. [8]) 
can be accommodated in this formalism. Since however we prefer not to rely on such 
specific assumptions, we will consider A (^) as a free parameter. 

To NLO, the (resummed) expansion of 7 obtained from Eq. ([I8|) is related to the 
unresummed expansion obtained from Eq. (|T^ by 



7(iV, as) = 7s (t) + «.7.. (t) + • • • , (22) 

where 

^« ("/v) ~ ^■^ (n-ax) ' 

Iss {jf) = Iss {n-A\) ~ ( ~as ^ V (^"^^ 

Since the resummation only involves formally subleading terms, 

ls + aslss = ls + aslss + 0{al/N). (24) 

A resummed double-leading expansion can finally be constructed by combining the 
resummed anomalous dimension 7 (|2^) with the standard expansion of 7 at fixed N. This 
gives a resummed double-leading expression for the large anomalous dimension eigenvec- 
tor. It can further be shown [13] that resummed double-leading expressions for the full 
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matrix of anomalous dimensions and for coefficient functions can be obtained by perform- 
ing the replacement N ^ N — AX in all remaining quantities, i.e. the projectors and the 
coefficient functions. 

The construction of the resummed double-leading expansion entails a further ambi- 



guity in the treatment of the double counting subtractions in Eq. (|16]): because these 
terms are common to the fixed-iV expansion 70, 71, ... and the fixed Ug/N expansion 
7s, 7ss, ... , we are free to decide whether to leave them unaffected by the replacement 
N N — A A or not. This defines a pair of resummation procedures, which of course 
only differ by subleading terms. Clearly, a variety of intermediate alternatives would 
also be possible. The main difference between these prescriptions is the nature of the 
small singularities of the anomalous dimension, which control the asymptotic small 
X behaviour. The resummed anomalous dimension always has a cut starting at = A 
Eq. (0), which corresponds [13] to an x^'^ behaviour of splitting functions at small x. If 
the subtractions are affected by the replacement, then 70 and 71 are the same as in the 
unresummed case (S-resummation), i.e. they have a simple pole at = 0, which leads 
to a "double-scaling" [19] rise at small x. If the subtractions are unaffected, this pole is 
removed by the subtraction itself (R-resummation). It follows that if A is positive, then 
the two resummations give similar results at small x, namely an x^'^ power rise. If A < 0, 
the S-resummation will display double scaling at small x, while the R-resummation will 
display a valence-like x"^ behaviour. 



4 Predictions for THERA 

A comparison of the resummation discussed in the previous sections to recent HERA 
data [20] was presented in ref. [13]. The best-fit results of that reference can be used 
to obtain predictions for THERA and discuss the study of small x scaling violations at 
such a facility. Because of the larger center-of-mass energy available at THERA, these 
predictions essentially amount to an extension of the current kinematic range of 1/x by 
about a decade for each value of Q^. It is interesting to note that more or less the 
same center-of-mass energy would be available at a hypothetical lepton-hadron collider 
obtained combining LEP with the LHC 

The phenomenological analysis of ref. [13] is based on a fit to data for the reduced 
cross-section 

aredix,y,Q^) = F2(x,Q2) + ^Fi(x,Q2). (25) 

2(1 - 2/) + r 

determined from structure functions F2 and computed to next-to-leading order in 
the double leading expansion with the R- and S-resummation prescriptions discussed 
in Sect. |^, by evolving parton distributions given at a scale Qq = 2 GeV. A standard 
unresummed next-to-leading order fit is also performed for comparison. The fits are 
performed in the parton scheme, so the quark distribution coincides by construction with 
F2. The large-x shape of parton distributions is taken from a global fit, while the small x 
behaviour is parametrized by two free parameters Xg and A^,, which give the asymptotic 
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Figure 2: The structure function F2{x,Q'^) obtained from a fit [13] to HERA data [20]. 
The prediction for THERA is the last decade in x for each value of Q^. The solid curve is 
an unresummed fixed-order two loop fit, while the dot-dashed curve corresponds to the 
S-resummation and the dashed curve to the R-resummation discussed in Sect. ^. 

small-x behaviour of the singlet quark and gluon distributions respectively as x^^^ . 
The resummation parameter A Eq. (|2T|) is also left as a free parameter. The strong 
coupling is fixed at ^^(M^) = 0.119. 

Because at the initial scale Qq abundant data are available down to the smallest values 
of X, and F2 coincides with the quark distribution, the quark exponent Ag turns out to be 
the same in all fits, and gives the effective power rise of the F2{x, Qq) ~ x~^'': Ag ~ 0.2. 

The best-fit value of the gluon exponent is valence-like in all fits: in the two-loop fit it 
is Ag ~ —0.1; while in the resummed fits it is significantly more valence-like, A^ ~ —0.2 
for both S- and R-resummation. The value of the resummation parameter A instead 
varies significantly according to the resummation prescription which is adopted. For the 
S-resummation, any value A < gives a good fit, with the best fit around A ~ —0.25. 
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Figure 3: Same as Fig. |], but for the structure function Fl{x,Q'^). 



As discussed in Sect. |^ the S-resummation with vanishing or negative A is closest to 
the unresummed fixed-order result. With the R-resummation, instead, only a fine-tuned 
value of A ~ 0.2 gives a good fit. This value of A turns out to be the same which one gets by 
fine-tuning the resummed anomalous dimension so that it be closest to the unresummed 
one in the HERA kinematic region [7]. Both resummed fits give a similar ~ 52 with 
93 degrees of freedom, to be compared to the unresummed value = 60. 

The structure function F2{x,Q'^) obtained in these fits is displayed in Fig. ||. It is 
apparent that, given the high precision of the HERA data, all curves, which give good fits 
to the data, are constrained to lie essentially on top of each other throughout the HERA 
region, except possibly at the smallest x values x ~ 10~^ at the initial scale Qo, where 
the R-resummation curve rises slightly less. In the THERA range it is still very difficult 
to tell the difference between various prescriptions at higher scales > 100 GeV^, but 
at lower scales, while results in the S-resummation are still essentially indistinguishable 
from the two-loop ones, the R-resummation predicts a somewhat faster evolution. 
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X 

Figure 4: Same as Fig. ^, but for the MS gluon distribution G{x, Q"^) = xg{x, Q^). 

The structure function Fl is displayed in Fig. This structure function is not deter- 
mined very accurately by the HERA data for the reduced cross section Eq. (^), essentially 
because of the scarcity of large-?/ data. The spread of the results is accordingly larger. 
Because at small x has a large gluonic component, the behaviour of Fl is similar to 
that of the gluon distribution, displayed in Fig. ^ For ease of comparison with other 
work, the MS gluon is shown, even though our fits were perfomed in the DIS scheme. 
Both resummations give a rather softer behaviour than the fixed-order one at the initial 
scale Qo- valence-like for the gluon, turned into a rise of Fl at very small x (well into the 
THERA range) by the rise of the coefficient function. The R-resummation, however, then 
leads to significantly more rapid evolution: as the scale increases, the resummed gluon 
overtakes the fixed-order one. This is essentially due to the fact that the R-resummation 
eventually generates an x~^ power behaviour of all parton distributions at small x, and 
here A = 0.2 (power rise). The S-resummation, instead, leads to evolution dominated by 
double scaling, which is very similar to the fixed-order one, and thus both the gluon and 
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Fl preserve the relative softness that they displayed at the initial scale. 

Summarizing, it is clear that resummation effects, though very small, become increas- 
ingly important as x decreases. Deviations from the fixed-order behaviour appear, at 
least in a simultaneous determination of F2 and F^. At present, the deviations from 
the fixed order prediction are within the uncertainties of the resummation procedure: so, 
while it is clear that the resummed gluon distribution is softer than the unresummed 
one, it is hard to tell whether it will evolve faster or slower at small x. The underlying 
physics between these options is quite different: either the onset of a slow power- like 
rise (R- resummation) , or persistence of the double-scaling rise (S-resummation). Both 
possibilities are consistent with present-day data, as well as with our current knowledge 
of anomalous dimensions. Understanding which (if any) of these possibilities is correct 
could be of considerable theoretical interest, and in particular, it could shed light on the 
running of the coupling in the high-energy limit [13, 16]. 

In conclusion, accurate data in the THERA region could reveal significant differences 
between the resummations procedures, and thus shed light on the structure of unknown 
higher order contributions to perturbative anomalous dimensions, and on the underlying 
physics. The simultaneous measurement of F2 and Fl in a wide range of Q'^ at small x 
would allow an accurate determination of structure functions at small x which are required 
e.g. for precise phenomenology of heavy quark production at future colliders. 

Acknowledgements: This work was supported in part by EU TMR contract FMRX- 
CT98-0194 (DG 12 - MIHT). 
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